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Abstract. One of the methods used to extend two-dimensional bosonization to 
four space-time dimensions involves a transformation to new spatial variables so 
that only one of them appears kinematically. The problem is then reduced to 
an Abelian version of two-dimensional bosonization with extra "internal" coor- 
dinates. On a formal level, putting these internal coordinates on a finite lattice 
seems to provide a well-defined prescription for calculating correlation functions. 
However, in the infinite-lattice or continuum limits, certain difficulties appear 
that require very delicate specification of all of the many limiting procedures 
involved in the construction. 

Introduction 

After the initial success of bosonization in 1+1 dimensions [1,2] Alan Luther 
presented a heuristic formula [3] for bosonization of free, massless relativistic 
fcrmions in 3+1 dimensions. A configuration-space transform [4] was later used 
by H. Aratyn [5] to put this into a prettier but no less heuristic form. For other 
approaches to this problem the reader is referred to References [6] and [7]. 

The basic idea of the method described here is to view field theories in 3+1 
dimensions as 1+1 dimensional theories with an internal, non-kinetic degree of 
freedom. The main subject of this paper is to examine carefully the limiting pro- 
cedures involved in bosonizing a 1+1 dimensional free fermion with an internal 
degree of freedom that takes values that are in a continuum. 

Tomographic transform 

We begin by reviewing the transformation to a 1+1 dimensional theory both for 
bosons and for fermions. We then present the standard Abelian technique for 
bosonizing the 1+1 dimensional fermion when the internal degrees of freedom 
lie on a finite lattice. There follows a careful analysis of the continuum limit of 
this lattice in the context of a simplified model. 



* Based on a talk given by C. Sommerfield at Giirsey Memorial Conference I, 
Bogazi§i University, Istanbul, Turkey, June 1994. 



The tomographic transform of a free scalar field </>(x, t) with mass to in 3+1 
dimensions is defined as 



4>(y, n, t) = y d 3 x - n • x)</>(x, i). 

Here n is a unit vector and cPn an element of solid angle in the direction of n. 

The field equation satisfied by <f>(y, n, t) and its equal time commutation 
relations are essentially those of a 1+1 dimensional scalar field, also free and 
with mass to. 

The corresponding transform for a fermion field is 



°(y,n,t) = ± J d 3 xd,%-n.x)^ U t a (n)^(x,t). 



It satisfies the 1+1 dimensional Dirac equation for a right-moving fermion field, 
(that is, depending only on y — t) and has the expected equal-time anticommu- 
tation relations. Here w a (n), a = 1,2, are orthonormal four-component spinors 
that satisfy (a ■ n)u a (n) = u a (n) where a are the Dirac alpha matrices. 

The strategy is to use standard Abclian 1+1 dimensional bosonization to 
construct the right-moving ip a (y, n, t), in terms of the right-moving part of the 
transformed boson (f)(y,n,t), assuming, to begin with, that n is restricted to a 
lattice. This is given formally, with the t dependence suppressed, as 

/oo 
dy' e(y-y')d t 4> a (y',n)}, 
-OO 

where e(y — y') = (y — y r )/\y — y'\. It follows that 

(0|C(y, n)$(y', n')|0) = - ^^(n, n') In (p[ a - i(y - y')]) 

where [i and a are infrared and ultraviolet cutoffs, respectively. 

We first take n to be discretely distributed on some lattice of directions and 
reinterpret the Dirac delta function as a Kronecker delta function. Replacing 
the pair a and n by a single index A, we find the formal bosonization expression 

where 

K A - exp \^Y. C eAC [^(°°) - <Pr(-™)\ 

and where e AB = —e BA , with le' 45 ! 2 = 1. The Klein factors K A are needed 
to make sure that fermion fields with different indices anticommute rather than 
commute. In what follows we will ignore these factors. Their only effect is to 
correct an occasional sign to its proper value. 



The correlation functions for the fermion field are determined from the 
matrix elements of their bosonic representation in the bosonic vacuum. The 
two-point function is computed as 



(0\i> A (y)ft B (y>)\0) = j 



2ir {n[a + i(y-y>)]y AB - 

If A = B we get the known fermion function. If A ^ B then this vanishes 
in the limit fi — > 0, so that the result is indeed proportional to S AB . When 
reinterpreted in terms of a continuous n we do get the correct two-point function 
for the transformed fermion field: 

(0|fr(y,n)ffV, n ')|0) = j- fT^'L 

2tt a — i[y — y'j 

All of the other fermion correlation functions come out properly as well. 

To test the consistency of the bosonization formula we must consider the 
important fermion bilinear operators in the 3+1 dimensional theory such as 
the Poincarc group generators and the chiral charge operators. If is one of 
these, does the sequence fi[^>3+i] <-> f2[-0i+i] <-> f2[</>i+i] <-» Q[<^3 + i] make sense? 
Trouble arises when we consider the rotation and boost generators. The proper 
treatment of the continuum limit of the lattice plays a fundamental role here. 

Simplified model and smoothing functions 

In order to get at the heart of the problem of the continuum limit we will treat 
a simple case in which the the internal variable is one-dimensional. We thus 
consider a single component 1+1 dimensional massless chiral fermion field ip(x, u) 
that depends on a continuous internal variable u whose domain is the real line. 
The bosonization of the fermion field in the Dirac equation (d t + d x )ip(x, u) = 
will be in terms of a right-moving chiral massless boson field (j> r (x, u). Assuming 
for the moment that u is a discrete variable, we write, up to Klein factors, 
ip(x,u) = (27ra) -1 / 2 exp [iy/^T</> r (x, u)] and obtain for the fermion two-point 
function 

m(x,u)^(x',u')\0) = ±- ^p-— 

Zir a — i(x — x') 

For the bosonization procedure to work, it is clear that the exponent that 
appears in the evaluation of the fermion two-point function must be associated 
with a Kronecker delta function, rather than a Dirac delta function so that 
the logarithmic singularity in the bosonic two-point function exponentiates to 
a simple pole or a simple zero. We need a more analytic method of converting 
one type delta function to the other. Returning to a continuous u, we introduce 
a new boson field <&{x,u) — J du'f(u— u')<fi r (x,u'). We take f(u) to be a real, 
even function. The Dirac delta function is replaced in the commutation relations 
among the $ operators, as well as in the bosonic two-point function, by the 
convolution g(u — u') = J du" f(u — u")f(u" — u 1 ). It is thereby softened. We 



will choose f(u) so that g(u) > 0, g(0) = 1, and e(u)g'(u) < if g(u) ^ 0. 
To make sure that g(u) goes to very quickly as |u| increases we take it to 
be a function of u/X and eventually take the scale parameter A to zero. This 
simulates the properties of a Kronecker delta function. 

We define the field ijj(x,u) = C(27ra) -1 / 2 exp [i\/47r$(x, u)] which will be 
the candidate for the canonical fermion field ip in the appropriate limit. Here C 
is some constant, to be specified later, that depends on /i and on the choice of 
g(u). We then find for the fermionic two-point function 



{OMx, utyV. «')|o> - \c\ 2 vP-°^± 



2ir [a - i(x - x')]9( u - u ') ' 

As the infra-red cutoff \i goes to zero, this vanishes unless u = u'. Choosing C to 
depend on \x so that C diverges suitably in this limit, we obtain the desired Dirac 
delta function and the same answer as in the discrete case, with the Kronecker 
delta replaced by the Dirac delta. 

There is a surprise when we consider the four-point function. One must 
analyze the behavior of the quantity |C| 4 // (2+fll2+fl34_Sl3_fl23_Sl4_fl24) as \i -> 0. 
Here gij = g(Ui — Uj) with u\ and U2 assigned to the fermion fields and and 
ii4 to the adjoint fields. In order to obtain the correct behavior when all four 
internal indices are close we find that we cannot choose f(u) such that g'(0) = 0. 
In fact g(u) must have a cusp at u = 0. Thus Gaussian functions for f(u) and 
g(u) are ruled out. A triangular pulse is acceptable for g(u) and this corresponds 
to a rectangular pulse for f(u). 

Fermion bilinears 

We go on to discuss the fermion bilinears. This is where we had trouble in the 
discrete case. A generic form of such a bilinear operator is [?/>' (x2, 112), ip(xi, Ui)] 
where the arguments are allowed to come together, possibly after various deriva- 
tives have been taken. This is the nature of the spatially point-split structure of 
the charge operator, and of the generators of translations in x, t and u. We pro- 
ceed to evaluate it using operator-product expansion methods (which are exact 
in this model). We have 

W2M = -^fi (1 - gi2) (R 12 - 912 - i? 2 r 921 ):e tv ^ ( * 1 "* 2) :. 

Here R\ 2 = a — i{x\ — x 2 ). The normal ordering is with respect to the bosonic 
creation and annihilation operators. We introduce u = \{u\ + u 2 ), v = m — u 2 , 
x = \{xi + x 2 ) and £ = x\ — x 2 , and obtain 

where fi(£,t>) = : cxp iy/4n[®(x + ^£,u + \v) - - \t,,u - \v)]i. Taking 
v and £ small, (while keeping a <C £), and sending (i to we discover that 

[$M->-i8(vMZ,0)/(nQ. 



The chiral charge operator Q is given in terms of the canonical fermions by 



Q = \j dx du\ du 2 5(u\ — u%)[^{x, 112), ip(x,ui)]. 

We choose to smooth out the u delta function and replace it by N\g(ui — u 2 ) 
where N\ is chosen so that N\ J du\ g(ui — u 2 ) = 1. Then N\g(ui—U2) becomes 
a Dirac delta function as A — * 0. We construct the candidate chiral charge Q in 
terms of ip as Q = lim^o N\ / dx du dv g{v)\\^}\^ V'l] so that 

Q = lim ~^ (0) Jdxdu [fi(f,0)-fi(-£,0)] = ^= J dx du ^{x,u). 

The charge density in u space is thus Q(u) — Nxtt^ 1 / 2 [$(00, 0) — $(— 00, 0)]. 
We then have [4>(x, u), Q(u')] — N\g(u — u')ip(x,u) which becomes the correct 
canonical relation as A — > 0. 

The canonical operator J = —i\ J dx du [tp^(x, u), jf^ip(x, u)] + (h. c.) gen- 
erates translations in u. It is the analogue of the rotation and boost operators in 
the tomographic representation of the 3+1 dimensional case. If we were to follow 
the procedure used for the chiral charge we would find that when expressed in 
terms of bosonic fields, the candidate generator would vanish. The remedy for 
this is to back up a step and delay taking u to zero. Then the spreading function 
used to define the bilinear, (which need not be the same as what we used for the 
chiral charge, namely N\g(u)), can be taken to depend on u in just such a way 
as to yield the correct answer in terms of canonical boson fields in the limit. 
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